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We discuss the Josephson current in superconductor / dirty normal conductor / superconductor 
junctions, where the superconductors have d^2_y2 pairing symmetry. The low-temperature behavior 
of the Josephson current depends on the orientation angle between the crystalline axis and the normal 
of the junction interface. We show that the ensemble-averaged Josephson current vanishes when 
the orientation angle is 7r/4 and the normal conductor is in the diffusive transport regime. The 
rfa.2_„2-wave pairing symmetry is responsible for this fact. 
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I. INTRODUCTION 



In recent years, the Josephson effect between the 
anisotropic superconductors has attracted much atten- 
tion because the high-Tf, superconductors might have the 
dx2-y2 pairing symmetryoD. In anisotropic superconduc- 
tors, the sign of the pair potential depends on the di- 
rection of a quasiparticle's motion. As a consequence, 
the zero-energy states (ZES)tl are formed at the normal 
metal / superconductor (NS) interface when the poten- 
tial barrier at the interface is large enough. The ZES 
are sensitive to the orientation angle between the crys- 
talline axis of the high- Tc superconductors and the nor- 
mal of the junction interface. The zero-bias conductance- 
peak which is due to the ZES is observed in conductance 
spectra of N / I / c?-wave superconductor junctions and 
the peak-beight is maximum when the orientation an- 
gle is 7r/4Q. Here I is the insulator. In SIS junctions, 
the ZES dominate the dc Josephson effect and the low- 
temperature anomaly in the Josephson currwit-has been 
discussed in a number of theoretical workg3il3. In ex- 
periments, it seems to be difficult to fabricate clean SIS 
and SNS junctions. Thus it is important to understand 
the effects of disorder on the Josephson current. So far it 
has been pointed out that the roughness at the interface 
suppresses the low-temperaturg anomaly of the Joseph- 
son current in SIS junctionsfflla. 

In this paper, we study the dc Josephson effect in su- 
perconductor / normal conductor / superconductor junc- 
tions, where the superconductors have the dj;2_j^2-wave 
pairing symmetry and the normal metal is in the diffu- 
sive transport regime. We analytically derive the Josepii^ 
son current in dirty SNS junctions based on a formulate 
in which the Josephson current is-|Calculated from the 
two Andreev reflection coefficientala. We show that the 
ensemble-averaged Josephson current vanishes when the 
orientation angle is 7r/4. The analytic results are con- 
firmed by a numerical spiulation by using the recursive 
Green function methodtJ. Throughout this paper we 
take units of fc^ = h = 1, where fc^ is the Boltzmann 



constant. 

This paper is organized as follows. In Sec. 2, we derive 
the general expression of the Josephson current in SNS 
junctions. In Sec. 3, we discuss the Josephson current in 
dirty SNS junctions. The discussion is given in Sec. 4. 
We summarize this paper in Sec. 5. 



II. JOSEPHSON CURRENT FORMULA 

Let us consider the two-dimensional SNS junction as 
shown in Fig. 0, where the length of the normal metal is 
ijv and the width of the junction is W, respectively. The 
pair potential in the momentum space is schematically 
depicted in the superconductors. In the dx2_y2-wave su- 
perconductor, the a-axis points the direction in which 
the amplitude of the pair potential takes its maximum. 
The angle between the x direction and the a-axis is aL 
{an) in the superconductor on the left (right) hand side. 
We assume that the effective mass (m) of an electron 
and the Fermi energy (/^f = kp/2'm) are common in the 
superconductors and the normal conductor, where fc^ is 
the Fermi wavenumber. We describe the SNS junction 
by using the Bogoliubov-de Gennes (BdG) equatiorll3. 
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where the first term in Eq. (^) is the barrier potential at 
the NS interface and the second term is the impurity po- 
tential in the normal metal, respectively. The energy of 
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a quasiparticle {E) is measured from the chemical poten- 
tial of the junction. The pair potential can be described 
by the Fourier representation, 



A(R, v-v') = j^j Jdk A(k) e'''^^ 



(4) 



where we assume that the pair potential is uniform in 
the superconductors and neglect the dependence of A on 
R — (r + r')/2. The s-wave superconductor is character- 
ized by A(k) — Ae"^J . The dx2-y2-wave superconductor 
is characterized by 



A(k) = Ae'^^- cos(2aj - 27), 



=,17 



57 + isin7 = [kx + iky)/kF, 



(5) 
(6) 



where A is the uniform amplitude of the pair potential, 
ipj and aj with (j = L or R) are the phase oipthe pair 
potential and the orientation angle, respectively^. At the 
Fermi surface, the wavenumber in the x direction is kx 
and that in the y direction is fcy, respectively. In the 
normal conductor, the pair potential is set to be zero. 

We calculate the JoseDhaon current by using the 
Furusaki-Tsukada formulaQ'E^, 



J 
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with i = L or R. Here T and u;„ = (2n -|- 1)ttT are 
the temperature and the Matsubara frequency, respec- 
tively. The wavenumber in the y direction ky specifies 
the propagating channel and the number of propagation 
channels for each spin is Nc = Wkp/n. The coefficients 
ai{ky,ujn) and a2(ky,ujn) are the analytic continuation 
{E —t iun) of the reflection coefflcients of a quasiparti- 
cle from the left superconductor to the left superconduc- 
tor. The Andreev reflection coefficient from the electron 
(hole) channels to the hole (electron) channels is denoted 
by ai (02)- In the presence of the time reversal symme- 
try, the Josephson current can be decomposed into the 
series of 



J : 



E 



(10) 



where ip = ip^ — ^Pr- In this paper, we neglect the series 
of Jm for TO > 2. This approximation is justified in the 
presence of the high potential barrier at the NS interface. 
In Fig. ^ (a), we show the four reflection processes of oi 
and 02 which contribute to Ji. In order to estimate oi 
and 02, we calculate the transmission and the reflection 
coefficients of the single NS interface for fixed ky. The 



sixteen coefficients are obtained from the continuity con- 
dition of the wavefunction at the NS interface since there 
are four incoming and four outgoing channels for each 
ky. In Appendix A, we show eight transmission and four 
reflection coefficients which are required in the following 
calculation. The Andreev reflection coefficients in Fig. Q 
(a) are estimated as 



al\ky,bJri) 



^t'^%{ky,aL,(pL) tk^,k'^ 



X r^Nik'y, an, (Pr) tt ,^^ t''j^s{ky,aL,^L), (11) 



af-'(fcj,,a;„) 



^ts%{ky,aL,ipL) tl k> 



X r%\{k'y, aR, <pr) t'^,^ ,^^ t%^s{ky,aL, <Pl), (12) 
4^^(fcy,w„) = Vi) ife„.fc; 

ai\ky,UJn) = ^tf^(ky,aL,LpL) y 



(13) 



X ''W7v(fcy, OIR, (Pr) tl,^ k^ t%s{ky, aL^ipL), (14) 

where tl^i'^l is the transmission coefficient of the elec- 
tronlike (holelike) quasiparticle in the normal conductor, 
and ky indicates the propagating channel at the right 
NS interface. The transmission coefficients in the normal 
metal are described by 
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where Qi^^ (r, r') is the Green function of the normal con- 
ductor, vp is the Fermi velocity, and Xky{y) is the wave- 
function in the y directipn belonging to the channel spec- 
ified by fcj,, respectivelyt3. In this paper, we assume that 
the NS interface is sufficiently clean so that ky is con- 
served while the transmission and the reflection at the 
interface. In aj^' in Eq. (pr]), a quasiparticle wave is ini- 
tially incident into the normal part from the left super- 
conductor through the channel specifled by ky . After the 
Andreev reflection at the right NS interface, we assume 
that the reflected wave transmits to the left supercon- 
ductor through the initial channel of ky because of the 
retro property of a quasiparticle under the time reversal 
symmetryll3. The two Andreev reflection coefficients in 



Eq.(^ are given by ai 



.(1) 



,(2) 



and a2 



^2 ^ + '^2 ^ J 



respectively. By using the transmission and the reflection 
coefficients in Appendix A, the following equations can 
be derived, 
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where (• • •) means the ensemble average and is the av- 
erage over the propagation channels. The transmission 
coefficients are approximately given by 



m JO 
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dy dy'X{LN,y;0,y'), (28) 



with 



X(r,r')^(t/.„(r,r')f;-.„(r,r')). 



(29) 



^tk,k,^ (19) 



(20) 



The Josephson current in Eq. (^ can be written as 

J = 2ieT^ ^ (P1-P2). (21) 

1^,1 ky.k'y 

The equation in Eq. (|21) corresponds to the reflection 
processes shown in Fig. | (b). After small algebra, we 
find the final expression of the Josephson current, 

J = 4esinv5T^ E Q^K^'^L)Q{k'y,aR)T{ky,k'y), 

l^-n ky,k'^ 

(22) 



In the diffusive regime, the function X{r, r') for small ujn 
satisfies the diffusion equation. 

To {2\co„\ - DoVl) X(r, r') « 27riVoTo(5(r - r'), (30) 

where Dq^ tq and Nq are the diffusion constant, the mean 
free time and the density of states at the Fermi energy 
per unit area for each spin degree of freedom, respec- 
tively. The propagating process in the normal conductor 
is diagrammatically illustrated in Fig. |^, where the cross 
represents the impurity scattering. We solve Eq. ( pO| ) and 
show the results in Appendix B. The averaged transmis- 
sion coefficients are 



In 



9N 



sinh/n 



(31) 



with 



where In = \J2n + lLfq/£^jj, = ^Dq/2ttT is the co- 
herence length, and Gat = {2e^ /h)gN is the conductance 
in the normal metal, respectively. The ensemble averaged 
Josephson current in dirty SNS junctions is calculated by 
using Eqs. (||) and ^ 



cos2 7A"i4:+ 
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{J{aL,aR)) ^ eTy^ N{aL)N{aR)gt, . , , , 
^ — ' smhm 
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(32) 



(33) 



— ^ (Aj Aj + ) + cos jA^ , (25) where we replace the summation J2k ^^"^ Integra- 



(26) 



where Z — mVb/kp represents the strength of the barrier 
potential. In what follows, we consider the high potential 
barrier at the NS interface, (i.e., Z ^ 1). We note that 
Q{ky,aj) is proportional to the Andreev refiection coef- 
ficient at the NS interface. Equations ( |2l] ) and (|2^) are 
the general expressions of the Josephson current propor- 
tional to simp. It is possible to apply these expressions 
to various junctions if the transmission coefficients in the 
normal part can be calculated. 



III. DIRTY SNS JUNCTIONS 



When the normal conductor is in the diffusive trans- 
port regime, i^/''], is almost independent of the propa- 



tion (Nc/2) J^^J^^^d-f cos 'J. When the orientation angle 
are zero, (i.e., a^^R = 0), the ensemble-averaged critical 
current in units of tt Aq / eRj is given by. 
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(J(0,0))--— T , 
25 Aq ^-^ ui + A-^ smhm 



where the resistance of the junction is described by 

-1 



(34) 



R.i = 



h QZ^ 



(35) 



gating channels. 



{T{kyX)) ^ {tH^), 



(27) 



and Aq is the amplitude of the pair potential at T = 0. 
In the case of the s-wave junctions, the numerical factor 
9/25 in Eq. (^) is replaced by 1 and the expression is 
identicai|ta that by the Usadel equatiorO. In the recent 
studiesEjlEil and references therein, it is pointed out that 
the results in Rcf. ^ is not correct in the low tempera- 
ture regime and the non-sinusoidal current-phase relation 
is observed. In these studies, the electron transmission 
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at the NS interface is perfect and the higher harmonics 
contribute to the Josephson current. In our results, the 
ampHtude of the Josephson current takes its maximum 
at if = 7r/2 for aU temperature regime because of the 
potential barrier at the NS interface. 

When aL,R = 7r/4, we find that iV(7r/4) = 0, therefore, 



(J(7r/4,V4)) =0. 



(36) 



This is because and '^l,r are even function of 

7, whereas — —A sin 27 is odd function of 7 in 

Eq. (I33I). The symmetry of the pair potential is re- 
sponsible for the disappearance of the ensemble-averaged 
Josephson current. 

When the orientation angle is tt/8, the averaged criti- 
cal current results in 



(J(7r/8,V8)) 



97r2 A 
512 A^" 



T 



In 



A2/2 sinh/n 



(37) 



The condition for the ZES (i.e., A+A^ < 0) is satisfied 
in the range of 7r/8 < I7I < 37r/8 in Eq. (|3^). In the same 
reason with the case of aL,B. = 7r/4, the contribution from 
this range becomes zero. In Eq. ( p7|) , the integration in 
the range of < I7I < tt/S in Eq. (33) contributes to the 
averaged Josephson current. We conclude that the ZES 
do not contribute to the averaged Josephson current in 
dirty SNS junctions of the d-wave superconductor. 

In Fig. U (a) and (b), we show the numerical results 
of the Josephson current at ip = Tr/2 by using the re- 
cursive Green function method on the two-dimensional 
tight binding modcllj, where Aq = 0.01/iF, = lOOao, 
W — 30ao and qq is the lattice constant. The mean 
free path {£) is about 6.2ao. The numerical results for 
"L,i? — and 7r/4 are shown in Fig. ^ (a) and (b), re- 
spectively. The dependence of A on T is described by 
the BCS theory. The lines are the Josephson current for 
several samples with different random configurations and 
the open circles denote the ensemble average. In (a), the 
results for all samples increase with decreasing the tem- 
perature and the averaged Josephson current agrees with 
the analytic results in Eq. (jsj) which are plotted in Fig. ^ 
(c), where a parameter Ln /£^r)\T=T^ = 7.6 is estimated 
from the numerical simulation. In Fig. ^(b), the averaged 
Josephson current is almost zero for aL^R = 7r/4 as it is 
predicted in Eq. (|3^). This fact, however, does not mean 
the absence of the Josephson current in experiments for a 
single sample. The sign of the Josephson current is either 
positive or negative depending on the random configura- 
tion and |J| of a single sample increases rapidly with 
decreasing the temperature as shown in Fig. ^(b) . These 
results indicate an importance oL-the mesoscopic fluc- 
tuations in the Josephson currentt3uj. In Eq. (p7|), we 
assume that {T{ky,k'y)) is independent of ky and k'y. In 
other words, the transmission coefficients are isotropic in 
the momentum space. In order to explain the numerical 
results in Fig. ^(b) , we have to consider a sample-specific 



anisotropy in the transmission coefficients. Here we in- 
troduce the anisotropy in a simple function. 



T{ky,k'y)^{tH'^) 

1 + /l Sky.ki 4^,fei + /2 Sky.k2 Sk'y,-k2 



(38) 



where fi and /2 are positive numerical factors much 
smaller than unity. In Eq. (^8|), we consider the situ- 
ation where the two elements T{ki,ki) and T(fc2,— ^2) 
are slightly larger than the average. By using Eq. (|3^), 
the Josephson current in a single sample with aL.B. = 7r/4 
results in 



J(,/4,,/4).|.Ti:-'" 



sinh/n 



X [fiYM - /2r(72)] , 

3 A cos^ 7 sin^ 27 

4 [|a;„| + cos2 7|sin27|]^ 



(39) 
(40) 



The denominator of Eq. (|o|) approaches to zero in the 
limit of Wn ^ and Z 3> 1, which reflects the ZES at the 
NS interface. In Fig. || (d), we plot the analytical results 
in Eq. ( |39| ) for several choices of /i , /2 , 71 and 72 , where 
/i and /2 are of the order of 10^^. Some of the results 
show non-monotonic temperature dependence, and oth- 
ers change the sign with decreasing the temperature as 
shown in both Figs. U(b) and (d). The analytical results 
well explain the numerical results. These results indicate 
the enhancement of the mesoscopic fluctuations of the 
Josephson current in the limit of T — > 0. 



IV. DISCUSSION 

In Sec. Ill, we calculate T{ky,k'y) by solving the diffu- 
sion equation in Eq. (|^. The conclusion in Eq. (^) is 
independent of the detail of the approximation. In the 
diffusive conductors, {'T{ky,ky)) are independent of ky 
and k'y. The conclusion in Eq. ( ^ ) always holds, when 
{Tiky, k'y)) is isotropic in the momentum space because 
the integration with respect to 7 at the two NS interfaces 
can be carried out independently in Eq. (p2[). For this 
reason, we derive an equation for the dirty SNS junctions, 

(J(«L,^/4)) - (J(^/4,a;^)> =0 (41) 

for any ul and ur. We also derive relations, 

{J{aL,aR)) = (J(q!l, -a_R,)) 
= ( J(-aL, afl)) = ( J(-aL, -aR.)) (42) 

because the zero energy states do not contribute to the 
ensemble average of the Josephson current. 

In the s-wave SNS junctions, the amplitude of the fluc- 
tuations is dJc — {eEc/h)Cd, where Ec — TiDq/L?^ is the 
Thouless energy of the normai, conductor and Cd is a 
constant of the order of unityE3. The amplitude of the 
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fluctuations in the Josephson current for aL,R = at 
T = are expected to be proportional to the Thouless 
energy and depend on the sample size since the charac- 
teristic features of the Josephson current are essentially 
the same with those in the s-wave junctions. On the 
other hand, the fluctuations for aL,B. — 7r/4 in Fig. ^(b) 
are mainly due to the ZES at the interfaces and strongly 
depend on Z and T as shown in Eq. ( ^0| ) and Fig. ^(d). 
Though the average of the Josephson current is zero, the 
fluctuations for aL,R = 7r/4 can be larger than those for 
C'L,R = 0. Further theoretical investigations are neces- 
sary to understand the amplitude of the fluctuations. 

The impurities in the normal metal may suppress the 
pair potentials near the NS interfaces because the nor- 
mal impurities break a Cooper pair in the anisotropic su- 
perconductors. The pair potential should be determined 
self-consistently in sweh situation. So far in SIS junc- 
tions, Tanaka, et. al.Ba calculated the Josephson current 
where the pair potentials determined in a self-consistent 
way and compared with the Josephson current obtained 
in a non-self-consistent manner. Their results show that 
there are no qualitative differences in the Josephson cur- 
rent between the two methods. In their study, the pair- 
breaker is the insulator. The pair-breaker in the present 
paper is the impurities in the normal metal. Though 
the origin of the pair-break are different in the two sys- 
tems, the suppression of the pair potentials, therefore, 
the suppression of the Josephson current is considered 
to be the common consequence when we determine the 
pair potential self-consistently. Thus it may be possible 
to infer that the Josephson current in real SNS junctions 
would be smaller than that of the present paper. The 
self-consistent study has to be done to discuss the am- 
plitude of the Josephson current quantitatively. However 
the main conclusion in Eq. ( p6| ) is not sensitive to the 
profile of the pair potential because the d-wave symme- 
try is responsible for the disappearance of the averaged 
Josephson current. The impurities may also modify the 
symmetry of the pair potentials. If the finite values of 
averaged Josephson current is observed in experiments, 
such results might be an evidence for the change of the 
pairing symmetry due to the impurities. This is because 
the d-wave component of the pair potentials does not 
contribute to the averaged current for a = 7r/4. 

The sign change of the pair potential on the Fermi sur- 
face, which is a characteristic feature of the anisotropic 
superconductors, leads to the disappearance of the 
Josephson current. Thus the conclusion in Eq. ( ^6| ) may 
be applied to the SNS junctions of the superconductors 
with p-wave pairing symmetry such as UPts for certain 
orientation angles between the junction interface normal 
and the node lines of the pair potentials. The investiga- 
tion in this direction is now in progress. 



V. CONCLUSION 

In conclusion, wc analytically derive the general ex- 
pression of the Josephson current in SNS junctions of the 
d-wave superconductors. In dirty SNS junctions, we show 
that the ZES do not contributes to the ensemble-averaged 
Josephson current because of the symmetry of the pair 
potential. In particular, when the orientation angle of 
the c?-wave superconductor is a = 7r/4, the ensemble- 
averaged Josephson current vanishes. The critical cur- 
rent of a single sample, however, remains finite value, 
which indicate an importance of the mesoscopic fluctua- 
tions of the Josephson current. 
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APPENDIX A: TRANSMISSION AND 
REFLECTION COEFFICIENTS 

We derive the transmission and reflection coefficients 
of the NS interface at a; = as shown in Fig. ||. In what 
follows, we calculate the coefficients after analytic con- 
tinuation (E — !■ iujn) for Un > and we assume that 
A <^ lip. In the normal metal, the wavefunction of the 
quasiparticle with ky = kp sin 7 can be described by 



■f3 



A 



B 



(Al) 



where a and /3 (A and B ) are the amplitudes of the in- 
coming (outgoing) waves in the electron and hole chan- 
nels, respectively. The wavefunction in the y direction is 
Xkyiv), V± = /cFC0S7±iw„/wFC0S7, and vp = kp/m is 
the Fermi velocity, respectively. 

In the same way, the wavefunction in the superconduc- 
tor is given by 



7$ 



xxfe„ (y), 



(A2) 



where 7 and 5 (C and D ) are the amplitudes of incom- 
ing (outgoing) waves in the electron and hole channels 
respectively. We define that 

A± 



A COS 2(7 q= aj), 



+ /2, 



(A3) 
(A4) 
(A5) 
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vf = /2, 

k± — kp cos 7 + ifl^ /vp cos 7, 



(A6) 
(A7) 



where j = L ov R. The phase and the sign of the pair 
potential is considered in the matrix, 






-i(¥>j+0±)/2 



(A8) 



where e"^-' = S(?n(A*). The two wavefunctions satisfy 
the continuity condition at the left NS interface (i.e., 
X — Q and aj ~ oll), 



dx 



(A9) 



(AlO) 



From Eqs. ( [A9| ) and ( |A10| ), the amplitudes of the out- 
going waves are connected with those of the incoming 
waves. 



(All) 



We explicitly show 10 coefficients which are required for 
calculating the Josephson current. 
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These coefficients are a function of fcy, and (p^. On 
the derivation, we use the approximation k± ~ p± ~ 
kp C0S7, where we assume that lu„ <ti fJ-p. 

The reflection coefficients at the right NS interface are 
obtained in the same way. The two Andreev coefficients 
are given by 



rNN{ky,aR, ipn) = -iQ{ky, aR)e 
r'}^N{ky,aR,ifiR) = -\Q{-ky,aR)e"^'^. 



(A29) 
(A30) 



APPENDIX B: SOLUTION OF THE DIFFUSION 
EQUATION 

In this paper, we consider the low temperature regime 
T < Tc ~ 10"^ - lQ~'^^p and the diffusive transport 
regime I/^ptq ^ 10^^ — 10^^. Thus arelation 27rTTo < 1 
is satisfied. In these regime, we solve the diffusion equa- 
tion 

TO {2\cjn\ - DoVl) X{r,r') « 2^Aoro(5(r - r'), (Bl) 

under the appropriate boundary condition0@. The 
right hand side of Eq. (Bl) corresponds to the first term 
in Fig. ^ which can be calculated to be 



^0 ^ 
vpr 



-rjl 



(B2) 



where I = vpTQ is the mean free path and r 



respectively. In Eq. (Bl), we replace this function by the 
(5-function in two-dimension because the smallest length 
scale in the diffusion equation is £. The function X{r, r') 
is the Green function of the Schrodinger-type equation 



DoV'fiv) = A/(r) 



(B3) 



The normal conductor is separated by the high poten- 
tial barrier at the two NS interfaces and there is no cur- 
rent flow in the y direction. In such situation, we solve 
Eq. (|B3|) under the boundary conditions, 



6 



dm 



dx 
df{r) 



= 0, 



dy 



x—O.Ln 



(B4) 
(B5) 



The function X(r,r') is represented by 

oo oo 

X{r,r')^2nNoY, E 

m=Q m'=0 
{"Tf) <^OS (^) COS (^^) 



2\LUn\ + Do[{mn/LNy + {m'Tr/Wy] 



for m = 
, , for m 7^ 0, 

form^O 
^, form 7^ 0. 



(B6) 
(B7) 

(B8) 



Next we calculate the function 



/ dy / dy' X{r,r') 
Jo Jo 



x—Ln,x'—0 



2ttNoTo 



W 1 



In 



Lff 2|a;„|ro sinhln 



(B9) 



where In = \J2n + IL^ /S^d and = \/Dq/2ttT is the 
coherence length. Finally the averaged transmission co- 
efficients are 



In 



1 



9N 



sinhZn 2|ct;„|ro ' 



(BIO) 



where — 2hN[)DqW / Lj^ is the dimensionless conduc- 
tance of the normal metal for each spin degree of free- 



dom. In Eq. ( BIO ), {f^t^) seems to be singular in the 
limit oi ujn 0. We note that the singula rity is stem- 
ming from the boundary condition in Eq. (B4). When 
there is no potential barrier at the NS interface, we ap- 
ply the boundary condition /(r)U=o,LN — instead of 
Eq. (B4) and obtain 



Ir 



9N 



sinh In 



3h(2^2|L^„|To) - 1 



(Bll) 



where the integration with respect to y and y' are carried 
out dX X = Lfq + y/2 and x' = —i/y/2. Since the prop- 
agation of a quasiparticle in the normal conductor is not 
sensitive to the boundary condition at the NS interface 
and N^{t'^t'^) is clo se to qn in the limit of w„ 0, we 



regularize the Eq. (BIO) by introducing the cut-off, 



Ir 



9N- 



sinh In 

F{x) =e(l-x) + e(a;- 



^^(2|c^„|to), 



(B12) 
(B13) 



where 8(a;) is the step function. Since {t'^t'^) decays as 
exp(— V2n -l- ILtv/^d), the contributions from small n 
are dominant in the summation with respect to Un in 
Eq. (E^). Thus in most cases, F(2|a;„|To) is unity. 
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FIG. 1. The SNS junction of the d-wave superconductor 
is illustrated. The orientation angle in the left (right) super- 
conductor is az, («/{). 

FIG. 2. The reflection coefficients in (a) contribute to the 
Joscphson current proportional to sm{ipL —tpR). The Joseph- 
son current calculated from the four reflection coefficients in 
(a) is summarized in a reflection process in (b). 



FIG. 3. The propagation process in the diff'usive normal 
conductor is diagrammatically illustrated. The cross repre- 
sents the impurity scattering. 

FIG. 4. The Josephson current in dirty SNS junctions is 
shown as a function of the temperatures. The numerical re- 
sults for aL.R = and 7r/4 arc plotted in (a) and (b), respec- 
tively. The lines are the numerical results of several samples 
with different random configurations and the open circles are 
the average over a number of samples. For comparison, we 
show the analytical results for aL,R = in (c) and for 7r/4 in 
(d), respectively. 

FIG. 5. The transmission and the reflection coefficients at 
the left NS interface. There are four incoming and outgoing 
channels. 
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